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THE PROBLEM 



When the particles of a fluid are endowed with some scalar density S, the 
density evolves in time according to 

DS 

dfS + u • VS = — - = thermal noise. (1) 

The right-hand side represents the microscopic spreading of S on the molec- 
ular level, and can be thought of as noise added onto the fluid velocity u. It 
is normally described by a term like KV'~.b' where k is a diffusivity. The study 
of (1), especially for chaotic flows and turbulent flows, has been extensively 
carried on for many decades. ^'^ 

Fluid motions also move vector flelds around. An evolving vector fleld V 
is governed by an equation of the form 

VY 

dtV + u • VV — V • Vu = = thermal noise. (2) 

The advective derivative in (1) is replaced by a Lie derivative VY /Vt in 
the evolution operator. The extra term enters because the evolution of a 
vector fleld involves not only being carried about by a fluid particle but 
also being turned and stretched by the motion of neighboring particles. The 
thermal noise on the right side is usually a simple diffusion term, at least for 
the case of constant density. Density variations bring in some unappetizing 
complications that we shall ignore. 

If the right sides of (1) and (2) are zero (perfect fluid motion) then S 
and V are frozen-in properties, and the fluid motions can distort any V in a 
complex fashion. In particular, when the dynamical system 

X = u (3) 

produces chaotic motion, the effect on the distribution of the advective fields 
can be rich and surprising, giving rise to intense local concentrations and 
lacunae of fine structure. 

In real fluid settings the situation may be more complicated than even 
this description suggests. If either 5* or V can feed back on the dynamics of 
u itself, the equations lose even their superficially linear appearances. For 
ordinary fiuid motions, we have pV = w, where p is the fiuid density and the 
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problem is no longer kinematic. Rather, u satisfies 

p{dtU +ujxu) — — Vp — ^pVu^ + F + thermal viscous effects, (4) 

where F is an external force density. So when pV is vorticity, we have coupled 
equations for u and V. 

G.l. Taylor, who made early contributions to the study of (1) for turbulent 
u, observed that vorticity is concentrated by turbulence.^ To learn what 
properties of the motion favor this effect, we may begin with the study of (2) 
without worrying, at first, whether these motions correspond to solutions of 
(4). That leads to the search for a field u that may produce local vorticity 
enhancement when introduced into (2). We call this the kinematic turbulence 
problem, after the usage of dynamo theory, where (2) applies when pV is the 
magnetic field. 

Batchelor and others have sought analogies between the vorticity and 
magnetic field, since they are both controlled by equation (2). However, this 
viewpoint has been belittled because u = Vxu, while no analogous relation 
exists between B and u. Yet this relation is implied by (2) and (4), and it 
need not as a result be considered exphcitly. When B is the field in question, 
wc couple it to (4) through a nonlinear term representing the Lorentz force. 
That is, 

p{dtu + u • Vu) = - Vp - -^B X V X B + F + thermal viscous effects. (5) 

47r 

Typically the kinematic aspects of the dynamo problem are stressed over 
those of the analogous vorticity problem. But if the analogy is good, it is not 
clear whether this should raise interest in the kinematic turbulence problem 
or throw the kinematic dynamo problem in a bad light. In any case, there 
does seem to be considerable interest in studying the effect of fiuid motion 
on an immersed vector field through (2), for prescribed u. 

These general thoughts lead into our discussion of the effect of a chaotic 
motion on a vector field satisfying (2). We describe in general terms the 
procedures that we have been following, omit the most strenuous portions 
of the calculations and refer to the related work on the mixing of a passive 
scalar by a continuous flow^, and to the study by Aurell and Gilbert^ of fast 
dynamos in discrete maps. 
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THE FORMAL SOLUTION 

We consider velocity fields u that are steady. In this autonomous case, 
the problem can be simplified by writing V(x, t) = Vo(x) exp(At). Then we 
obtain, for the diffusionless case, 

VV/Vt = (u • Wo - Vo-Vu)e^* = -AVoe^*. (6) 

We may look for a solution of the form Vq = g(x)u(x), where 

u • Vg = -Xq. (7) 

This shows how the kinematic problem for vector fields may be related to 

the more extensively studied problem of passive scalar transport^; if we set 
S'(x, t) = g(x) exp(At) in (1), we get (7). Moreover, if q is constant on 
streamlines, then we must have A = 0. These special solutions arise in the 
very restricted conditions of steady flow without dispersive effects, and they 
illustrate the kind of degeneracy that we encounter when the conditions are 
too simple. 

More generally, we can write the solution of (2) formally, as shown by 
Cauchy. Let x(t, a) be the position of the fluid particle that was at the point 
a when t = 0. Then the field evolves according to 

V(x,i) = J(a,i)V(a,0) , (8) 

where J'(a, t) = d{x.)/d{a) is the Jacobian matrix of the transformation that 
moves the fluid into itself with time with x = x(a, t). 

We write x = y>*a, where <p^ is the flow that maps the initial positions of 
the fluid particles into their positions at time t. Its inverse, a = i^~*x, maps 
particles at time t and position x back to their initial positions. Then we 
can write (8) in the seemingly complicated, but quite useful form, 

V(x, t)^ J S{a - v-*x) J(a, t) V(a, 0) d^a. , (9) 

where the integral operator introduced here is the analogue of the Pcrron- 
Frobenius operator for the case of scalar advection.^ Having turned the dif- 
ferential equation (2) into an integral equation, we may use analogues of 
Fredholm's methods^ to solve it. 
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If we were to include the effects of diffusion by using a term of the form 
r/V^V on the right hand side of (2), with r] constant, we would need to replace 
T>/T>t by V/Vt — 77 in (6) and modify (9), but for now we shall speak only 
of the noise-free case. 

STRUCTURES OF CHAOTIC FLOWS 

To describe the evolution of a frozcn-in magnetic field, we need a suitable 
means to characterize the flow. For the chaotic flows on which we concen- 
trate here, the Lagrangian orbits determined by (3) densely contain unstable 
periodic orbits. Their union forms what we may call a chaotic set. If object 
is attracting, the term strange attractor is appropriate while in the opposite 
case it becomes a strange reseller, from which representative points depart 
as time proceeds. Even though all periodic orbits within such an ensemble 
are unstable, trajectories of the flow often spend extended periods tracking 
these paths. The strange set provides a delicate, skeletal structure embedded 
within the flow that can be used to systematically approximate such prop- 
erties as its local stability and geometry.^ This is the basic idea underlying 
cycle expansions^, the method that we apply here to the problem of vector 
advection. 

To be explicit, we consider a trajectory x{t) generated by the ordinary 
differential equation 

x + x-cx + x^ = , (10) 

with parameter c. This is a special case of equations arising in multiply 
diffusive convection^. The absence of a second derivative in (10) ensures 
that the flow is solenoidal, in line with the kinds of flow most commonly 
studied in dynamo theory, in spite of the fact that every putative dynamo is 
compressible. For certain values of c, homoclinic orbits exist and we let c be 
close to such a value and, in particular, one such that in the neighborhood of 
the origin the flow satisfles Shil'nikov's criterion for the existence of an inflnity 
of unstable periodic orbits. Then, an inflnite sequence of intertwined saddle- 
node and period-doubling bifurcations creates a dense chaotic set, and the 
motion of points in the flow is chaotic. 

Details of the structure of this flow are presented elsewhere; here we 
mention only that x{t), under the conditions mentioned, is a sequence of 
pulses. Moreover, the separation of the k-th and {k — l)-st pulses, A^, may 
be expressed in terms of the two previous spacings. This provides us with a 
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timing map, which is a two-dimensional map resembhng in form the Henon 
map^^. Thus, to good approximation, 

Ak^f + an , (11) 

where T is a mean period, a is a small parameter, and an irregular timing 
fluctuation satisfying, 

Tk+i ^ I - arl - Tk-i , (12) 

which is the orientation and area preserving form of the Henon map. 

The form of any particular pulse is quite close to that of the homoclinic 
orbit, which can be computed at the outset where the map (12) determines 
a sequence of pulse positions. Once these are known, we can generate a 
complete and reasonably accurate solution for the velocity field. Thus we can 
reconstruct the entire flow from this simple map. Moreover, this map also 
contains the invariant information about the periodic orbits of the chaotic 
set (cycle topology and stability eigenvalues) and serves as a powerful tool 
in the construction of cycle expansions. 

In the following sections, we describe the technique of cycle expansions for 
the problem at hand. This needs no explicit specification of the velocity field 
u, but in the concluding section we report some numerical results obtained 
for the particular flow modeled by (12). 

EVOLUTION AND TRANSFER OPERATORS 

A frozen-in vector field is stretched, squeezed and swept around by the 
chaotic flow, evolving as described by (8). To compute the large-scale evolu- 
tion of the held, we rewrite (9) as 

y,(x,t)= /(i3a4(x,a)y,(a,0) , (13) 

with a kernel 

4(x,a) = 5(a-v-*x)|^ , (14) 

where summation over repeated indices is understood. The kernel con- 
trols the evolution of the embedded vector field. This transfer operator is 
linear and possesses (for nice hyperbolic systems, the so-called Axiom A 
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flows) a sequence of eigenvalues, e"'^"*, e~'^^^, e~'"^*, .... For large times, the ef- 
fect of £* is dominated by its leading eigenvalue, e~'^°* with -Re(z/o) < R^{^i)-i 
i = 1,2,3, .... In this way the transfer operator furnishes the fast dynamo 
rate, v = —pq. 

The operator Cjj was introduced in ref. 8 in order to study the "stability 
of strange sets" , and applied to discrete map models of fast dynamos in ref. 4. 
Here we apply it to continuous time flows. 

A TRACE FORMULA 

To calculate the leading eigenvalue of £*, we evaluate the trace 

tr(£*)=y^d3a4(a,a) = y^d3a5(a-v-*a)^ , (15) 

which asymptotes to e"^ for long times. We evaluate this integral by means 
of explicit periodic orbit expansions. Each cycle within the fabric of the 
flow contributes to the overall trace of the operator, and each contribution is 
obtained by integration along the periodic orbit in suitable local coordinates. 
Wc omit the details of this purely technical operation (see ref. 3) and report 
that each prime cycle p together with its repeats contributes a term 

oo ir 

T y E S(t - rTj,) (16) 

^tt|det(l-J-)| ^ ^ ' 

to the integral. Tp is the cycle period and 3p is the transverse stability matrix 
u(t + Tp) = Jpu(t) for a two- vector in the tangent plane transverse to the 
flow. The Jp eigenvalues Ap i, Ap 2 are independent of the position along the 
orbit and the choice of transverse coordinates. 

The trace of the transfer operator is the sum over all periodic orbit con- 
tributions, with each cycle weighted by its intrinsic stability 

00 j-j. Y 

^^^•)-ZT.Zj^i^^^^^i*-rT.) (17) 
where we have introduced the Fourier representation of Dirac delta functions. 
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A FREDHOLM DETERMINANT 

We rotate the axes by the Wick transformation, k is, and observe 
that the summations in (18) correspond to the logarithmic derivative of the 
function 

°° 1 tr j: 



F(s) = exp 



-EE 



p ^-^|det(l- J-O 
related to the above trace by 



■e 



srTr, 



(19) 



tr(£*) = - / ds e-^'F'{s)/F{s) . (20) 

The function F{s) is a Fredholm determinant. Its zeros produce singular- 
ities in the integrand of (20). If it is an entire function (defined over the whole 
complex plane), then the contour of integration can be suitably deformed so 
as to encircle the various poles of F'(s)/F{s), and the trace becomes a sum 
over the residues of the integrand: 

oo 

tr(£*) = E^ne-^"* ' (21) 

n=0 

where is a pole of multiplicity m„. Hence, the spectrum of the transfer 
operator is determined by the zeros of the Fredholm determinant F{s). 

In order to simplify F{s), we factor the denominator cycle stability de- 
terminants into products of expanding and contracting eigenvalues. The 
example at hand is a 3-dimensional hyperbolic flow with cycles possessing 
one expanding eigenvalue Ap (of absolute value > 1), and one contracting 
eigenvalue Ap, with |Ap| < 1. Then the determinant may be expanded as 
follows: 

-j^ oo oo 

|det(l-J7)r =|(l-A7)(l-A7)ri = |AprEEAr''Aj'- . (22) 

j=0 k=0 

With this decomposition we can rewrite the exponent in (19) as 

oo 1 M^' I \^'\p^^'^p °° oo 1 

T. T. - /I '\ = E i: i: ; VAje"-')' (a; + a;) , (23) 

p r=l ^ I dct I 1 - J M I p j^k=0 r=l ^ 



which has the form of the expansion of a logarithm: 

E E [log (1 - e^^l A,| AJ-^ Aj) + log (1 - e^^l A,| A;UJ+^)] . (24) 

p j,k 
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The Fredholm determinant is therefore of the form, 

= Fe(s)Fe(s) , (25) 



where 



with 



Fe(s) = n n (i - i^'^K) , (26) 

P j,k=0 

CO 

Fc{s) = l[ n (1-4^'^ A,), (27) 

P j,k=0 

tW^e^T,\X,\^ . (29) 



The two factors present in F{s) correspond to the two Floquet exponents 
of the periodic orbits, the expanding and contracting exponents. In form, 
they are the Selberg-typc products known in other areas of mathematical 
physics. ■'^^ 

DYNAMOS FOR MAPS VS. FLOWS 

For 2-d Hamiltonian volume preserving systems A = 1/A, and (26) re- 
duces to 

^e(^)=nn ' ■ (30) 

Denoting the eigenvalue sign by ap — Ap/|Ap|, the Hamiltonian zeta function 
(the j — k — part of the product (27)) is given by 

ycdynis) = n (i - ^p^'^') ■ (31) 

p 

This is a curious formula: the Hamiltonian zeta function depends only on the 
return times, not on the eigenvalues of the cycles. Furthermore, the identity 

A + 1/A 2 



|(1-A)(1-1/A)| |(1-A)(1-1/A)| 

substituted into (25) leads to a relation between the vector and scalar advec- 
tion Fredholm determinants: 

Fdynis) = Foisf/Cdynis) . (32) 
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The Frcdholm determinants in this equation are entire for nice hyperbohc 
(axiom A) systems, since botli of them correspond to multiphcative operators^ 
For maps with finite Markov partition the inverse zeta function (31) reduces 
to a polynomial since the time Tp = Up is an integer, and the curvature 
terms^ in the cycle expansion vanish. For example, for maps with complete 
binary partition, and with the fixed point stabilities of opposite signs, the 
cycle expansion reduces to 

l/Cdyn{s) = 1. (33) 

For such maps the dynamo Fredholm determinant is simply the square of 
the scalar advection Fredholm determinant, and therefore all its zeros are 
double. In other words, for maps, the fast dynamo rate equals to the scalar 

advection rate^^. 

However, for flows the dynamo effect is distinct from the scalar advection. 
For example, for flows with finite symbolic dynamics grammar, (32) implies 
that the dynamo zeta function is a ratio of two entire determinants: 

l/Qyn{s) = Fayn{s)/F^{s). (34) 

This relation implies that for flows the zeta function has double poles at the 
zeros of the scalar advection Fredholm determinant, with zeros of the dy- 
namo Fredholm determinant no longer coinciding with the zeros of the scalar 
advection Fredholm determinant; the leading zero of the dynamo Fredholm 
determinant is larger than the scalar advection rate, and the rate of decay of 
the magnetic field is no longer governed by the scalar advection. 

NUMERICAL INVESTIGATIONS 

For our numerical investigations of the kinematic dynamo effect we have 
used (11) with T = 1, a = 6 and a e [0, 1]. For a = 6 the Henon map has 
a repeller with complete binary grammar. We have computed all periodic 
orbits up to length 13. The period of the prime cycle p is obtained by 
summing up the contributions (11): 

np-l 

Tp^Up + a ^ Xp^i . (35) 

where Xp,i are the periodic points in the prime cycle p. In table 1 we list the 
stabilities and Xp^i for the cycles up to length < 6. From such listing one 
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can generate the prime periods for any parameter value a, feed this data into 
cycle expansions, and extract a set of leading eigenvalues by determining the 
zeros of F{s) in the complex plane by standard root finding routines such 
as the Newton-Raphson method. A visualization of the complex function 
F{s) is afforded by a contour plot of In |i^(s)|. In such plot the eigenvalues 
correspond to the minima of ln|F(s)|, and the range of validity of a cycle 
expansion is indicated by absence of fine structure beyond some Re(s). 

Figs. 1(a) and 1(b) show the contour plot of Predholm determinants 
for the scalar and vector advection respectively, for a — 0. For this case 
(constant return time) the dynamics is equivalent to a map, and the dynamo 
Predholm determinant should have double zeros coinciding with the zeros of 
the scalar advection Predholm determinant. In this computation we used all 
prime cycles up to topological length 12, with polynomial truncation yielding 
the 12 leading zeros in both cases. Since the zeros of the dynamo Predholm 
determinant are double zeros, in the polynomial truncation only the first few 
are accurate, with the nonleading terms converging more poorly. 

The result for maps should be contrasted to the result for flows: figs. 
2(a) and 2(b) show the corresponding Predholm determinants for a flow with 
large dispersion of return times, a — 1. Here the zeros of the two Predholm 
determinants do not coincide, and the leading zero of the dynamo Predholm 
determinant yields the dynamo rate. Parenthetically, this lifting of eigenvalue 
degeneracy is very noticeable already for return time dispersion as small as 
a — 0.01. Pig. 3 shows the contour plot corresponding to the dynamo zeta 
function (31). The leading zeros coincide with the leading zero of the dynamo 
Predholm determinant, but this zeta has a double pole at the escape rate zero 
of Fo{s), see eq. (34), which is very apparent in this plot. In contrast, the 
zeta function (33) for a discrete time map is trivial. 

Note that the dynamo zeta function does not require evaluation of cycle 
eigenvalues Ap, and even so the pole in the zeta yields the escape rate; we 
can compute the escape rate from the cycle periods Tp only by locating 
the first pole of the dynamo zeta function. This indicates an intimate and 
not yet elucidated connection between the periods and stabilities of general 
dynamical systems. 

Another interesting property of the fast dynamo flows is that in these 
systems the magnetic field can grow exponentially although the fiow is re- 
pelling, the escape rate is positive, and scalar advection density decreases 
exponentially. This is illustrated by the difference in the sign of the leading 
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eigenvalues in figs. 2(a) and 2(b). 

CONCLUSIONS 

We have introduced a new transfer operator for chaotic flows whose lead- 
ing eigenvalue yields the dynamo rate of the fast kinematic dynamo and 
applied cycle expansion of the Frcdholm determinant of the new operator to 
evaluation of its spectrum. The theory has been tested on a normal form 
model of the vector advecting dynamical flow. If the model is a simple map 
with constant time between two iterations, the dynamo rate is the same as 
the escape rate of scalar quantities. However, a spread in Poincare section 
return times lifts the degeneracy of the vector and scalar advection rates, 
and leads to dynamo rates that dominate over the scalar advection rates. 
For sufficiently large time spreads we have even found repellers for which 
the magnetic field grows exponentially, even though the scalar densities are 
decaying exponentially. 
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FIGURE CAPTIONS 

Figure 1: Contour plot of log \F{s) \ for (a) the scalar advection Fredholm 
determinant, and for (b) the dynamo Fredholm determinant for the Henon 
map a = 0, using prime cycles up to topological length 12. Note that the 
scalar and vector advection rates coincide. 

Figure 2: Contour plot of log |-F(s)| for (a) the scalar advection Fredholm 
determinant for the flow (35), and (b) the dynamo Fredholm determinant for 
the a = 1 flow, using prime cycles up to topological length 12. Note that 
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for flows the vector advection eigenvalue (fast dynamo rate) dominates the 
scalar advection eigenvalue. 

Figure 3: Same as fig. 2 for the dynamo zeta. The leading eigenvalue 
yields the fast dynamo rate and the convergence radius of the cycle ex- 
pansion is controlled by the pole at the scalar advection eigenvalue eq. (34). 
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Table 1. All periodic orbits up to 6 bounces for the Hamiltonian Henon 
mapping (11), a — 6. The columns hst the topological length of the cycle, 
its expanding eigenvalue Ap, the period of the orbit, and the binary code for 
the cycle. 
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